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A simplified analytical model of tropospheric random variations in angle 
measurements is described. This model is used to predict the minimum and 
maximum power density spectra between which the tropospheric random 
angle errors of observations on the Early Bird satellite are expected to lie. 

The apparent angular position of the Early Bird satellite was then meas- 
ured at microwave frequencies with the large horn-reflector antenna at the 
AT&T station near Andover, Maine. Random variations in the azimuth 
and elevation angles have been observed and recorded. The analysis of these 
records results in a description of the observed random angle variations by 
their power density spectra. 

A comparison of the predicted power density spectra from the model with 
the observed spectra is made. It is concluded that the observed random angle 
variations are indeed caused by random tropospheric refraction. 

The feasibility of acquiring data on atmospheric propagation effects, 
particularly tropospheric angle errors, with the aid of geostationary satellites 
is therefore also demonstrated. 

I. INTRODUCTION 

1.1 Objective of this Paper 

The performance of earth-based radar and optical systems is ulti- 
mately limited by temporal and spatial random variations in the refrac- 
tive index of the tropospheric propagation medium. It is the objective 
of this paper to present a method for predicting random tropospheric 
angle errors in such systems, and to compare a prediction with micro- 
wave observations made on the Early Bird satellite. 
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1.2 Problem Approach 

The scintillation or twinkling of the stars which is experienced in ob- 
servations through the earth's troposphere is a familiar effect of the 
random variations in the refractive index of this propagation medium. 
Astronomers have known for a long time that the troposphere actually 
causes variations in at least four characteristics of the received star light, 
namely: (i) the intensity, (it) the spectral distribution of the intensity, 
(iii) the shape of the telescopic diffraction image, and (iv) the apparent 
angular position of the star. Scientific studies 1 - 2 of these effects seem to 
concentrate mainly on the intensity scintillations. The random varia- 
tions in the other characteristics, especially in the apparent angular posi- 
tions of stars, are treated in much less detail. 

However, in those radar and optical systems which are used to measure 
the position (and its time derivatives) of both distant and near objects 
(such as aerospace vehicles) the random tropospheric angle variations 
assume great importance. For the analysis and synthesis of these sys- 
tems, it is valuable to accumulate the knowledge on the random tropo- 
spheric errors in form of a sufficiently general model. 

Such an analytical model of tropospheric random errors in the position 
measurements and then time derivatives in radar and optical systems 
has been developed. Among other capabilities this model also permits 
the prediction of the random tropospheric angle variations for specified 
sets of tracking situations and system parameters. The prediction is 
made in terms of minimum and maximum power density spectra (PDS) 
between which the observed spectra are expected to lie. 

The choice of PDS for the characterization of the random errors is 
necessary because the relation between errors at two points in this sys- 
tem is usually a function of the error frequency (/). The resulting PDS 
further permit (i) subsequent studies of the effects of frequency de- 
pendent data processing operations (smoothing, calculation of deriva- 
tives, prediction, etc.). (ii) detailed comparison with errors from other 
sources, and (iii) application of the optimization methods described by 
H. W. Bode, C. E. Shannon, and S. Darlington. 3 Values for the more 
familiar variance (a 2 ) or standard deviation (a) of these random errors 
at the output of these processes may then be obtained with a straight- 
forward integration of the output PDS (see Section 1.4 below). 

Within the model, the predicted PDS of the random tropospheric 
angle errors are analytically calculated by operating with certain model 
functions on a model power density spectrum winch is given in the range 
coordinate. This range model PDS represents the pooled data on tropo- 
spheric random refraction. It is based upon observations of random 
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variations in the tropospheric refractive index, and in range and phase 
measurements mainly made at the National Bureau of Standards. 4 • 5 - 6 - 7 

The successful launch of the Communications Satellite Corporation's 
Early Bird Satellite on April 6, 1965, and its subsequent stabilization in 
an almost perfect geo-stationary orbit, provided an opportunity to test 
the model. For this purpose, azimuth and elevation angle measurements 
on the microwave beacon of the Early Bird satellite were made with the 
large horn-reflector antenna at the American Telephone and Telegraph 
(AT &T) Station near Andover, Maine. Most of this ground equipment 
was previously described in detail. 91011 . 12 - 13 - 14 . 15 A brief description of the 
Early Bird satellite may be found in Ref . 16. 

The resulting angle error measurements are particularly valuable for 
comparison with the theoretical model since they are obtained under two 
unique conditions provided by a geo-stationary satellite as a target. 
First, the propagation path goes through the entire atmosphere so that 
the observed angle errors include possible effects of high altitude turbu- 
lence, which are impossible to obtain with Earth based targets. Second, 
the angular tracking rates are negligible relative to the effective wind in 
the troposphere with which the refractive index anomalies pass through 
the propagation path. 

Thus, the analysis of the angle error measurements and the prediction 
of the expected random angle errors for the geo-stationary satellite are 
considerably simplified compared to the analysis and prediction for the 
more frequent aerospace targets (aircraft, missiles, low satellites) which 
have large apparent angular velocities and motion disturbed by forces 
unknown in the necessary detail. The effects of temporal random varia- 
tions of the refractive index in the Earth's atmosphere on the angle 
measurements, integrated along the line-of-sight between a ground an- 
tenna and a geo-stationary satellite should be observable in an almost 
pure form. 

l .3 Scope of this Paper 

In this paper, a simplified version of the model of random tropospheric 
errors is first described, which permits the calculation of the predicted 
minimum and maximum PDS of the tropospheric angle errors for track- 
ing tasks involving one almost stationary point target and a single ob- 
server (single-site radar). 

Next, the presented model is used to calculate the numerical values of 
the predicted PDS of random tropospheric angle errors for the specific 
tracking situation of the Early Bird observations. 

The methods and specific circumstances of data acquisition for one 
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twenty minute period of observations of the Early Bird satellite from the 
AT &T ground station near Andover, Maine are then described. Another 
section is concerned with data processing and analysis; it includes the 
time series of observed azimuth and elevation angles, the calculations of 
their power density spectra and confidence limits, and estimates of the 
manual chart reading error and of the effects of thermal receiver noise. 

Finally, a comparison is made between the PDS of random tropo- 
spheric angle errors predicted with the model and the PDS of the ob- 
served random angle variations. 

1.4 Scaling of Power Density Spectra 

In this paper, the random variations of the observed azimuth and ele- 
vation angles will be described by their power density spectra (PDS). 
The numerical computation of the PDS from the time series of data is 
made by the indirect method described by Blackman and Tukey. 17 
However, the scaling of the PDS in this paper deviates from that of 
Blackman and Tukey by denning the variance <r 2 of the random error as 

* 2 = fp[f)df. (1) 

Jo 
Thus, the PDS P{f] is valid only for positive frequencies, / ^ 0. The 
power spectral density is 

P\f\ = ?$ (2) 

of the variance contribution d(<r 2 ) to the random error, per unit frequency 
bandwidth, df, at the frequency, /. 

II. PREDICTION OF RANDOM TROPOSPHERIC ANGLE ERRORS 

2.1 Model Concept 

The analytical model of random tropospheric errors in radar and op- 
tical systems which has been developed permits the calculation of the 
power density spectra, and variances of range, phase, range difference, 
and angle errors, and their time derivatives from a basic pool of model 
data with the aid of certain model functions. This general model ac- 
commodates many different sets of system parameters, and is flexible 
enough to allow modification for its continuing improvement based upon 
the analysis of additional data. 

During the development of the model the usual lack of sufficient data, 
and the non-stationarity of the tropospheric refractivity field soon made 
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themselves felt. It was realized therefore, that only an approximate 
model of the real troposphere could readily be constructed, which neces- 
sarily would yield approximate predictions. However, it was found that 
this approximate model was good enough to allow the useful prediction 
of tropospheric errors in several interesting cases of tracking system 
analysis and synthesis. 

In the following part of this paper a simplified version of the general 
model of random tropospheric errors is described, which is limited to the 
prediction of the power density spectra of the random tropospheric errors 
in the angle measurements made by a single-site radar (or radio tracker) 
on an (almost) stationary target. 

The particular coordinate of radar measurements selected for the col- 
lective description of the pool of basic model data is the slant-range 
coordinate. In this approach, all available observations of random tropo- 
spheric errors are first normalized to certain model conditions, and trans- 
formed into the slant-range coordinate. An analytical power density 
spectrum (PDS) at the lower limit of these normalized and transformed 
observations is then denned as the model PDS in range, P m {f) , where / 
is the (error-) frequency. 

The derivation of the PDS for the random tropospheric angle errors, 
and for other than the model conditions, is then achieved by processing 
the range model PDS, P m {f], with certain power gain functions, called 
the model functions. These model functions depend on such parameters 
of the tracking situation as the antenna diameter, slant-range and eleva- 
tion angle of the target, weather, and wind. 

It may be noted particularly, that in this simplified version of the 
model it is not necessary to do any explicit processing in the space do- 
main. Based upon the assumption of an isotropic, and frozen turbulence 
field of refractive index anomalies in the troposphere, all processing is 
confined to the frequency {/} -domain. 

The entire model also can be used in an inversion of the computational 
flow to yield, from new observations, additional information on the basic 
range model PDS, P m {f\, and on the model functions. 

2.2 Assumptions and Limitations 

The simplified analytical model of tropospheric random errors in radar 
and optical systems is subject to a number of assumptions and limita- 
tions : 

(i) The model is intended to yield tropospheric errors in tracking 
tasks where one point target is directly observed within the local horizon 
along a line-of -sight (LOS) by a single observer. 
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(ii) It is assumed that the random errors are small, thus the model 
functions are linear in the sense of being independent of the magnitude 
of the errors. 

This assumption is justified by the finding that the random errors in 
the observed quantities (range, angles) have relative magnitudes of only 
one part in 10 5 , or so. 

(Hi) It is assumed that the random errors are stationary during the 
calculation, or observation of one PDS. The spatial and temporal non- 
stationarities of the random tropospheric errors are only considered by 
the introduction of the "global" weather functions. Local anomalies, as 
well as diurnal and seasonal variations of the tropospheric random errors 
thus are not separated here. It is believed that more detailed knowledge 
in this respect is better obtained by direct measurements under the 
particular local circumstances of actual tracking situations. 

(iv) It is assumed that the random errors due to the tropospheric 
anomalies can be treated as if they were caused by the motion of a locally 
isotropic field of "frozen" turbulence through the line-of-sight with an 
effective wind speed (u) normal to the LOS. 

(v) The wavelength (X) of the transmitted electromagnetic waves is 
assumed to be small, say X ^ 10 [cm], in order to avoid basic theoretical 
difficulties which are manageable only if X <3C I, where I is the characteris- 
tic length of the tropospheric anomalies. 1 - 2 This assumption is also im- 
portant in order to avoid the effects of random propagation through the 
ionosphere. 

(vi) The size of the antenna system is small with respect to the diame- 
ter of the earth (flat earth assumption). 

(vii) The size of the antenna system is small enough to avoid the lack 
of correlation between the refractivity anomalies at large distances on 
the surface of the earth. 

2.3 Model Functions 

2.3.1 Model Power Density Spectrum in Range, P m {f) 

The conditions to which the available data 4 ' 8 ' 6 ' 7 on random tropo- 
spheric range, phase, and refractive index variations are normalized are: 
(i) effective tropospheric path length = L m = 15 [km]; 
(ii) effective wind speed normal to LOS = u m = 1 [m/sec]; 
(Hi) surface refractivity = N m = 10 6 (n m — 1 ) = 313, this value is 
the U.S. average, 8 and n m is the equivalent refractive index; 

(iv) known effects of variations of the surface refractivity are not 
corrected during data acquisition. 
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After transformation to the selected range coordinate the model power 
density spectrum, P m {f\, is then derived as an analytical approximation 
to the lower limit of all observations. 

The derived model PDS in range consists of five branches, which are 
linear in a log (power density) versus log (frequency) plot, namely: 

f9.6 X 10 +19 / +2 [m 2 /Hz] 

for 0g/g2.5X 10~ 8 [Hz] 



(3) 



1.5 X KTVtmVHz] 

for 2.5 X 10~ 8 ^ / ^ 1.0 X 10 5 [Hz] 

. ,4.7X10- U r 2 ' 5 [m 2 /Hz] 

1 ' for 1.0 X 10~ 5 £ / £ 1.0 X 10" 3 [Hz] 

1.5 X 10- 12 r 3 [m 2 /Hz] 

for 1.0 X 10~ 3 If^LOX 10 +2 [Hz] 

1.5 X 10- 6 r 6 [m 2 /Hz] 

for 1.0 X 10 +2 ^ / ^ «> [Hz] 

where the frequency / is to be inserted in hertz. This PDS is plotted 
in Fig. 1. 

2.3.2 Angle Scale Function, S a 

The PDS of random tropospheric angle (a) errors for a single antenna 
radar are obtained from the range model PDS by operating on P m {f] 
with the angle scale function, S„ . The derivation of this function is 
based upon the fact that refractivity anomalies of characteristic length 
I, or of wavenumber k, which drift through the LOS with the effective 
wind speed u,„ cause random error components of frequency 

/ = ujl = KU m /2ir. (4) 

To simplify the analysis the circular antenna aperture of diameter d 
is now approximated by an interferometer system of equal angle ac- 
curacy and baseline length 



B = 0.626 d 



(5) 



which lies in the plane of the angle being measured. The angle measure- 
ment is thought to be indirectly obtained by a range-difference (or phase- 
difference) measurement across the effective baseline length B. The 
tropospheric refractivity anomalies disturb this range-difference measure- 
ment to an amount that depends on the characteristic length I and on the 
antenna diameter d. 
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Fig. 1 — Model power density spectrum, P m {f] , of tropospheric random errors 
in the range coordinate versus error-frequency,/, in a log-log plot. Pm\f\ is valid 
for the model conditions in Section 2.3.1. The tropospheric anomalies have the 
characteristic length I. 

It is found that for relatively short characteristic lengths 

I ^ h = 2d (6) 



which cause the high error frequencies 

/ ^ /i = u m /k = 0.5 u m /d 



(7) 



the random range errors due to these anomalies at the two ends of the ef- 
fective baseline length B are practically uncorrected with each other. 
Thus, they cause a power density of the random error in the range-differ- 
ence measurement that is twice as large as that of the random error in a 
single range measurement. Analytically, this finding may be expressed 
with the aid of a range-difference scale function 

S* R = 2 for / ^ /i = 0.5 ujd. (8) 

The tropospheric refractivity anomalies with characteristic lengths 
larger than the critical length k , namely 



I ^ h = 2d 



(9) 
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cause the lower error frequencies 

/ ^ /i = u,„/h = 0.5 ujd. (10) 

In this frequency region, the induced random range errors at the two 
ends of the effective baseline B are more and more correlated as the char- 
acteristic length is increased. It is found that with respect to the range- 
difference errors across B the antenna behaves like a high pass filter with 
break frequency /i . Analytically, the resulting reduction in the power 
density of the low frequency random range-difference errors may be ex- 
pressed by another branch of the range-difference scale function, namely 

Sa r = 2(///!) 2 for / ^ /i = 0.5 ujd. (11) 

The multiplication of the range model PDS, P m {f\, with the two 
branches of Sar in their respective frequency regions would result in a 
PDS for the random tropospheric range-difference errors across the base- 
line B under model conditions. 7 

The last step in the derivation of the desired angle scale functions is 
based upon the assumption of small angular deviations relative to the 
axis of the antenna system. Then the angle error a is simply related to the 
range-difference error AR and the effective baseline length B by 

a = AR/B. (12) 

In terms of power densities, this relation permits the calculation of the 
angle scale function S a from the range-difference scale function Sar and 
B, in general, as 

S a = S LR /B\ (13) 

The combination of (5), (8), (11), and (13) finally yields the angle 
scale functions in two branches, namely 

S a = 20 (f/u m ) 2 for Og/S/il 

r (14) 

S a = hi tt for fi^f £ oo J 

The breakfrequency of the angle scale function is 

/i = 0.5 ujd, (15) 

where u m = 1 [m/sec] is the model wind speed taken normal to the LOS 
and in the plane of the angle measurement, and d is the diameter of the 
circular antenna aperture. 
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2.3.3 Aperture Smoothing Function, <£<, 

The spatial smoothing on tropospheric random error components that 
are due to refractivity anomalies of small characteristic length (Z <$C d) 
is another function of the antenna diameter, d. Here the combined effects 
of several small refractivity anomalies tend to cancel across the antenna 
aperture, hence the antenna acts like a low-pass filter in the error-fre- 
quency {/} -domain. A simplified aperture smoothing function which 
analytically represents this effect is 

* a = 1 for g / ^ / 2 1 

\> (16) 

*« = (M) for / 2 S/g coj 

where 

ft = 2 ujd (17) 

is the breakfrequency for aperture smoothing in angle measurements. 

2.3.4 Effective Path Length Function, A 

From theories of propagation through a uniformly turbulent random 
medium, as for example given by Chernov, 1 it is known that the vari- 
ance of phase, range and related errors is proportional to the path length. 
This proportionality holds in both the near-field and the far-field regions 
of the "scattering" refractivity anomalies of a given characteristic length, 
I. Therefore, it is possible to account for an effective tropospheric path 
length, L, which is different from the model path length L m , by multi- 
plying the power density with the effective path length function 

A = L/L m for g/ ^ oo. (18) 

The required effective tropospheric path length, L, may be calculated 
by integrating over the geometrical length differentials along the LOS, 
which are weighted with the square of the local average refractivity at 
the height of each layer of the atmosphere. The necessary data on the 
variation of the refractivity with height have been taken from Bean and 
Thayer. 8 If the height of the target is /i 2 ^ 10 [km] above the surface of 
the Earth, and the apparent elevation angle is E a ^ 3°, the effective 
tropospheric path length becomes 

L = L /sin E a , (19) 

where 

L = l 1 lL) (6.61 - 0.0LV s )[km] (20) 



ffij* 



TROPOSPIIERIC ANGLE ERRORS 



1449 



is the effective height of the troposphere for surface ref inactivities 

250 ^ N s ^ 450, 

and N m is the model surface refractivity. 

In this case, which is relevant to many radar and optical tracking 
tasks, one thus has the effective path length function as 

A = [ N °/ Nm l (6.61 - 0.01tf.)[km] for ^ / ^ ». (21) 
L m • sin hj a 

2.3.5 Weather Functions, W 

As stated in paragraph 2.3.1 the model power density spectrum, 
P m \f] , is denned as the lower limit of the available observations normal- 
ized to the model conditions. Essentially all (say 99 percent) normalized 
observations exhibit larger errors than given by P m . Consequently, all 
PDS directly derived from P m [f} for other coordinates and tracking 
situations also would only give the expected minimum errors. Since it is 
frequently desired to state more about the expected distribution of the 
derived PDS above the expected minimum level we have introduced cer- 
tain power gain functions, called weather functions, W g . These are de- 
fined as the maximum weather function TF mftx {/} which covers the maxi- 
mum errors previously observed, and the median weather function 
W me d\f}. On a "global" basis (actually only embracing all circumstances 
of previous observations entered into the model data), it is expected that 
50 percent of the measured PDS will lie above and below the PDS pre- 
dicted with TF me d {/}, and essentially all (say 99 percent) of the measured 
PDS will lie below the PDS predicted with W m „\f\ . 

The maximum weather function derived from available observations ' 



5,6,7 



IS 



W max \f\ = 



6 for ^ / ^ 1.00 X 10 5 [Hz] 

1.89 X 10 +8 / +1 ' 6 

for 1.00 X 10" 5 ^ / ^ 2.23 X 10~ 5 [Hz] 

20 for 2.23 X 10" 5 ^ / ^ 1.00 X 10~ 3 [Hz] 

6.32 X 10 +2 / + °' 5 

for 1.00 X 10 -3 ^ / ^ 1.00 X 10 _1 [Hz] 



200 for 1.00 X 10 ! Sf < oo 



where / is in hertz. 

The median weather function is taken as 



[Hz], 



(22) 
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^med{/) = (^ mnx {/})°- 5 . (23) 

Both functions are plotted in Fig. 2. Note that W m in = 1 by definition. 

2.3.6 Effective Wind Functions, U 

The purpose of the effective wind functions, U, is to introduce other 
magnitudes of effective wind speed, u a ?± u m , into the model. The deriva- 
tion of these effective wind functions rests upon the assumption that an 
isotropic, frozen turbulence field of refractivity anomalies exists in the 
troposphere which moves through the LOS with a constant effective 
wind speed component, u a , normal to the LOS and in the plane of the 
angle coordinate a. With this assumption, a given anomaly causes an 
angle error of a magnitude that is independent of u a , and of a frequency 

io 3 



r. 



-8 -7 -6 -5 -4 -3 -2 -t I 2 3 

LOG 10 (f/[Hz]) 

Fig. 2 — Maximum and median weather functions, Wmix{f] and Wmed{f\, 
respectively, versus decadic logarithm of the error frequency, /. 

that is proportional to u a . It was found that a PDS that is given in 
{f,u m } -space as a sum of branches of the form 

P a '{f,u m } = P '(f/f o y for / m . min ^ / ^ /„,, max (24) 

is transformed into an equivalent PDS in \f,u a ) -space by the relations 

P a "{f,u a ) = U P -P a 'if,u m \ 

for 
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(25) 



where the effective wind function for the transformation of the power 
density is 

Up = (u m /u a ) y+1 (26) 

and the effective wind function for the transformation of the frequency 
regions 
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U f = u a /u m . (27) 

In these relations the/ a , m i„ and/ a , ma x are the limits of the frequency 
region in which P a " is valid after the transformation to [f,u a ] -space. 

In the special case of small angular velocity of the LOS, the equivalent 
wind due to the angular rate is negligible compared to the natural winds 
in the atmosphere. The effective wind speed is then simply 

u„ = w a , (28) 

where w a is that component of the natural wind which is normal to the 
LOS and in the plane of the angle a. In this plane of the angle a, the 
atmospheric refractivity anomalies, on the average, appear to move 
through the LOS with the speed w a . 

The calculation of the effective wind speed for azimuth (a — » A) 
angle errors depends on the geometrical relations between the LOS and 
the natural average wind vector, w, Fig. 3. If A is the azimuth angle of 
the LOS, and 8 is the azimuth of the wind vector, their difference 

= 8 - A (29) 

can be used to calculate the effective wind speed for azimuth angle er- 
rors 

u A = \wa\ = | $0 | ■ | sin j8 |. (30) 

With the horizontal LOS component (see Fig. 3) 

ivi = \w\- cos/3 (31) 

the effective wind speed for elevation errors similarly becomes, Fig. 4, 

u s = \w B \ = | wx I • I sin E» | (32) 

or 

u E = | w | • | cos • sin E a \. (33 ) 

Only the magnitudes of the effective wind speeds are of interest in this 
special case of small angular velocity of the LOS, and a single-site radar. 

2.4 Computation of the Predicted PDS 

In the prediction of the PDS of random tropospheric angle errors for a 
particular tracking situation numerical values are inserted for all inde- 
pendent parameters in the model functions given above. The range model 
PDS, P m [f], is then multiplied by the model functions, within the limits 
of the stated frequency regions, in the following sequence: angle scale 
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NORTH 

TARGET 




OBSERVER' 

Fig. 3 — Horizontal projection of LOS, wind vector w, and azimuth angle A. 

function, aperture smoothing function, effective path length function, 
weather function, and effective wind functions. 

For the tracking situation of the Early Bird observations on May 7, 
1965 the numerical values for the parameters and model functions are 
given in the Appendix. The operation of these models functions on the 
range model PDS, P m \f\ , resulted in four PDS, namely a predicted mini- 
mum spectrum and a predicted maximum spectrum for each of the two 
angle coordinates azimuth and elevation. The resulting PDS P^.min , 
Pa, ma* and /Vmin , Pe, max are plotted over the interesting frequency 
range in Figs. 12, and 13. 

III. OBSERVATIONS OF RANDOM TROPOSPHERIC ANGLE ERRORS ON THE 
EARLY BIRD SATELLITE 



3.1 Data Acquisition 

3.1.1 Method and Equipment 

For the acquisition of the data on random angle variations in azimuth 
and elevation, the apparent position of the microwave beacon (frequency 
about 4000 MHz) of the Early Bird satellite was measured with the horn 




OBSERVER 



Fig. 4 — Projection of wind vector into vertical plane through LOS with 
apparent elevation angle E a ■ 
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antenna (aperture diameter d = 67.7 ft) and its associated equipment. 
During these measurements the communications carrier of the satellite 
was switched off; this resulted in an increase in beacon signal strength 
to such a level that the thermal receiver noise in the obtained angle 
measurements was negligible compared to the desired tropospheric 
random errors. 

The signal flow through the major pieces of equipment which were 
used is illustrated in Fig. 5. After acquisition of the Early Bird satellite 
beacon in the main beam (beamwidth d = 0.225 deg) of the horn an- 
tenna, (he antenna control was turned over to the vernier autotrack 
system, and the servo loop opened by switching off the hydraulic drive 
motors. The antenna was now fixed in an orientation indicated by the 
digital display of the azimuth (A) and elevation (E) angles given in 
degrees, and derived from digital data pickoff units, which have a pre- 
cision of encoding 11 of 0.00275[deg]. 

The satellite now appeared to drift through the fixed horn antenna 
beam in an irregular motion, which was partially due to motion in its 
true position (orbit), but also due to the refractive index variations in 
the intervening atmospheric propagation medium, and possibly other 
disturbances. The apparent angular position of the satellite relative to 
the electrical axis of the horn antenna on the ground was determined by 
the autotrack system which contains angle error sensing and processing 
equipment. The azimuth and elevation error signals, AA{t] and AE{t\ 
from the autotrack system were passed through low-pass recording filters 
before recording either by oscilloscope and camera, or by analog strip 
chart recorder. 

The photographic pictures of the oscilloscope display giving AA vs 
AE were only used for inspection. The strip chart recordings giving the 
AA \t\ , AE\t\ time series, however, were used for the more detailed analy- 
sis of the data, as described later. 

3.1.2 Propagation Path and Mean Satellite Motion 

During the measurements the propagation path pointed from the horn 
antenna near Andover, Maine, to the Early Bird satellite approximately 
at an azimuth angle A tt 128.5° (southeast), and an elevation angle 
E £& 24.5°. The slant range between ground antenna and satellite was 
about 24,300 [statute miles] ~ 39,100 [km]. The terrain surrounding the 
Earth Station may be described as a shallow bowl of perhaps 10-miles 
diameter surrounded by hills of up to about 3.5 [deg] elevation. 

The mean apparent satellite motion with respect to the azimuth and 
elevation angles given above consisted of (i) a small linear drift with an 
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Fig. 5 — Flow diagram of data acquisition. 

azimuth component of Ai = —2.08 X 10~ 3 [deg/hr], and an elevation 
component of E\ = —1.22 X 10~ 3 [deg/hr], plus (ii) a diurnal elliptical 
motion with peak-to-peak amplitudes of A = 0.266 [deg] in azimuth and 
$ = 0.245 [deg] in elevation. The net result of these components appears 
at the Earth Station as a slow motion of the satellite along a helical path 
seen under an oblique angle. This picture of the mean apparent satellite 
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motion was obtained by plotting the hourly azimuth and elevation angles 
from the digital display for a few days before the analyzed random angle 
error data were recorded. The random azimuth and elevation angle er- 
rors, AA and AE, which are subjects of this paper are superimposed on 
this mean apparent motion. 

3.1.3 Date and Time of Observations 

The random angle error data recorded on strip charts, and analyzed 
in this paper were taken on May 7, 1905 between about 23 hr:38 min 
EDT and 23 hr:59 min EDT. 

3.1.4 Weather Conditions 

Weather data at the horn antenna of the Andover Earth Station were 
not taken. However, the weather data may be estimated from those 
taken at a private station in nearby Rumford, Me. This estimation yields 
the following data: 18 Cloud cover 9/10, wind South 19 [statute miles/ 
hr], dry bulb temperature 41.8 [°F], dew point 35 [°F], and pressure 
28.5 [inches] = 965.0 [millibars]. 

3.1.5 Recording Filters 

The low-pass recording filters mentioned in Section 3.1.1 above were 
simple two-section RC filters. Since the source impedances feeding these 
filters are small, and the load impedances connected to their outputs are 
large compared to the resistances in the RC sections of the filters, their 
inverse power gain is 

F = a~ l = 1 + (a>T) J ((co7') 2 + 7). 

In this equation, F is the ratio of input power to output power, 
T = RC is the time constant of one filter section, and w = 2irf, where/ 
is the frequency. 

The power density of the random angle errors before the filters may 
then be obtained by multiplying the power density of the recorded 
random angle errors with the inverse power gain F. The filters which 
were used in these observations allowed a choice between two cutoff 
frequencies. The results of numerical calculations of the inverse power 
gains versus frequency for the "LOW", and "HIGH" filters are plotted 
in Fig. 6. 

3.1.6 Calibration 

The sensitivities of the recorded error voltages (after the recording 
filters) to errors in the azimuth and elevation angles with respect to the 
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Fig. 6 — Inverse power gain, F, of recording filter versus frequency, /. Fi 
for filter in "LOW" range. h\ for filter in "HIGH" range. 

electrical axis of the horn antenna were obtained by direct calibration 
on the Early Bird satellite. For this purpose, the antenna servo system 
was disabled, and manual angle offsets were then inserted and their 
effects on the strip chart records were measured. 

3.1.7 Oscilloscope Displays 

Photographs of oscilloscope displays of the random elevation error 
(A-E 1 ) versus the simultaneously occurring random azimuth error (AA) 
were also made. 

The photo tracing in Fig. 7 was obtained at 22:30 EDT May 7, 1965 
while the recording filters were in the "HIGH" range, and the exposure 
time was five seconds. It is obvious that in this sample of the higher fre- 
quency errors the peak-to-peak azimuth variations (about 40 micro- 
radians) are considerably larger than those of the elevation errors (20 
microradians). 

The photo tracing in Fig. 8 was taken at 22 : 38 EDT on the same date 
with the recording filters in the "LOW" range. The exposure time was 
two minutes. In this sample of the lower frequency errors the peak-to- 
peak azimuth variations (12 microradians) are slightly smaller than the 
elevation variations (18 microradians). 

As mentioned before, these photos were only used for inspection and 
not for numerical analysis. 
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Fig. 7 — Tracing of oscilloscope photograph of random elevation error (AE) 
versus random azimuth error (A.4). Time is the parameter. Recording filters in 
"HIGH" range. Exposure time: five seconds. 

3.2 Data Processing and Analysis 

3.2.1 General Methods and Equipment 

The data on azimuth (AA ) and elevation angles (AE) versus time were 
recorded on a strip chart recorder with the recording filters in the "LOW" 
range. The time series of azimuth and elevation angles were manually 
digitized at two second intervals. 

After the manual digitizing process the time series of azimuth and 
elevation variations were punched into cards for subsequent processing 
on the 7094 digital computer. 



3.2.2 Time Series of Observed Angle Variations 

The time series of the azimuth (AA) and elevation (AE) angle varia- 
tions are shown in Figs. 9 and 10, respectively. 

The total observation time was somewhat above twenty minutes. 
This observation time was limited by the mean apparent drift of the 
satellite in the fixed antenna beam. This drift resulted in the recording 
traces going off scale after a certain time. 

A total of 720 azimuth data points, and 666 elevation data points were 
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Fig. 8 — Tracing of oscilloscope photograph of random elevation error (AE) 
versus random azimuth error (AA). Time is the parameter. Recording filters in 
"LOW" range. Exposure time: two minutes. 

recorded. Due to the systematic drift the elevation record went off the 
recording scale sooner than the azimuth record. 

3.2.3 Power Density Spectra of Observed Angle Variations 

The random variations of the observed azimuth and elevation angles 
will also be described by their power density spectra (PDS) for compari- 
son with the predictions. The numerical computation of the PDS from 
the time series of data is made on a digital computer by the indirect 
method described by Blackmail and Tukey. 17 It proceeded in the follow- 
ing steps: calculation and removal of the mean, and of the linear trend 
in the series; tapering the first 5 percent (start) and the last 5 percent 
(end) of the time series with a cosine function; computation of the auto- 
correlation function versus number r of 2£ r ^ M time lags each of 
duration of the sampling period At; computation of the Fourier transform 
of the autocorrelation function by a cosine series resulting in a raw power 
spectrum and subsequent smoothing of the raw spectrum by sliding, 
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Fig. 9 — Time series of observed azimuth angles, AA, versus time, t. 

weighted averages of values for three neighboring frequency steps with 
weights 0.25, 0.50, and 0.25. 

The computer program actually calculates a quantity X'\f\, called 
"power spectrum", which is related to the usual power density spectrum 
P'\f\ by the equation 



where 



p'\n = x'if)/Af, 



Aj = f N /M = l/(2-At-M) 



(34) 



(35) 



and 

f N = Nyquist frequency 

M = maximum number of lags in autocorrelation 

At = sampling period. 

In this equation, the primed quantities indicate that they still refer to 
the data at the output side of the recording filter. In order to obtain the 
desired power density spectrum at the input of the recording filter, 
P'\f\ must be multiplied by the inverse power gain of the filter, F{f] , 
yielding 

P{f\ =2-At-M-F{j)-X'{j). (36) 

Additional smoothing of the power density spectrum is used at the 
higher error-frequencies, since many cycles of these angle error compo- 
nents have been observed. This is done with a filter of approximately 
constant relative bandwidth, /3 = b/f = 0.231, at the expense of absolute 
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Fig. 10 — Time series of observed elevation angles, AE, versus time, t. 

frequency resolution. This point will be illuminated again in Section 
3.2.4. 

The data on azimuth and elevation angle variations given in Section 
3.2.2 were analyzed with the methods just described. It was found that 
the mean linear trends during these observations were in azimuth 
+8.1 X 10~ 8 [rad/sec], and in elevation -1.7 X 10" 7 [rad/sec]. Even at 
a distance of 10 [km] along the line-of-sight the magnitude of these angu- 
lar rates amount to beam sweeping speeds of less than 0.002 [m/sec], 
which are indeed negligible compared to natural wind speeds in the 
troposphere. 

The power density spectra of the observed azimuth and elevation 
variations at the input of the recording niters, Pa{S\ and Pslf}, which 
result from these calculations are plotted in Figs. 12 and 13, respectively. 
Other spectra also plotted in these figures are explained below. 



3.2.4 Confidence Limits for Power Density Spectra 

In computing confidence limits for the power density spectra it is 
necessary to distinguish between two error-frequency regions: the low- 
frequency region in which the absolute analyzing bandwidth of the PDS 
calculation 



b = (At -At)' 



(37) 
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is constant, and the high-frequency region in which the relative analyz- 
ing bandwidth 

= Wf (38) 

is constant. 

In the low-frequency region, with constant absolute analyzing band- 
width, b, the number of degrees of freedom in the PDS estimate is ap- 
proximately 

k = 2N/M (39) 

where N is the number of data points observed. 

In the high-frequency region, with constant relative analyzing band- 
width, /3, the number of degrees of freedom is frequency dependent ac- 
cording to 

k = 20JN-AL (40) 

The confidence limits for the calculated PDS of the observations can 
now be given, the lower limit being 

P 1 {f\ = Pifl/KJfl (41) 

and the upper limit 

Pitfl =K 2 \f}-P{f\, (42) 

where P\f) is the calculated PDS, and Ki,z[f) are the confidence factors. 
For a confidence level of p = 95 percent the upper confidence factor 
is approximately 

the total confidence factor (here only used as an intermediate to obtain 
Ki) 

( 240 \ / N 

K n = antilogio 1 /, = J , (44) 

and the lower confidence factor 

K 1 = K 21 /K 2 , (45) 

where k is the number of degrees of freedom given above. For k ^ 5 the 
stated analytical approximations for the confidence factors have less 
than 10 percent error. 

The resulting numerical values for the upper (if 2 ) and lower (Ki) 
confidence factors are plotted in Fig. 11. The results of calculating the 
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Fig. 11 — Relative confidence factors, Ki and K 2 , for p = 95 percent con- 
fidence level versus frequency, /. 

upper confidence limits (P A 2 , Pbi), and the lower confidence limits 
(P A \ , Pm) for the azimuth (P A ) and elevation (P E ) spectra are plotted 
in Figs. 12 and 13. 

3.2.5 Chart Reading Error 

The errors which are introduced into the data by the manual reading 
of strip chart records (digitizing) are of the same type as quantization 
errors. The variance due to a given quantization step size (q) is known 
to be 19 



<r, 2 = <?7l2. 



(46) 



If it is now assumed that the quantization noise, which causes this 
variance, is sharply bandlimited white noise of constant power density 
(Pq ), and with a cutoff frequency equal to the folding frequency of the 
digitized time series (Jn), then one also has the variance as 

»/- f N Pa-df= P g '-f N . (47) 

Consequently, the noise power density due to the manual chart reading 



is 



P,' = °Hh = ? 2 -A«/6. (48) 

As before, this primed power density is taken at the output side of the 



TROPOSPHERIC ANGLE ERRORS 



1403 



io-w 



io-' 2 
e 

6 
4 



IO"' 3 
10 







■\ 




















V 


















P A2 














\ 














P A , MIN^> 


■~-\ 












\ 


















P^ 




\ 
\ 








Pa,max\ 
















Pa, 








\ ^ 


























X 


































\ 


' \ 


^. 


\ 








V 




















\ 


^X 1 


\ V 








\ 




















\ 


+/~ 




*"N 






\ 


\ 






















\ 
\ 


_ < 


\ VV 






\ 


























\ N V 


. 






























Cm; 
































n3 






\ 


























/ 


\ 






























/ 


\ 






























/ 


\ 






























r 


\ 
















Pth 


















\ 
































\ 






















PqA 












\ 






























\ 






























\ 






























\ 












































\ 



































FREQUENCY IN HERTZ 



Fig. 12 — Power density of random azimuth angle errors, P, versus frequency, 
/. P A , P A i , P A > = observed PD!S and 95 percent confidence limits; pA.mia , 
Pa, max = predicted tropospheric PDS limits; P, h = thermal receiver noise; 

/',,,, = manual chart reading error. 

low-pass recording filter. In order to obtain the power density spectrum 
of the chart reading error referred to the input of the recording filter, it 
is necessary to multiply P q ' with the inverse power gain Fi\f] of the filter, 
see Section 3.1.5, which yields here for the azimuth coordinate 



JUf\ - HS\-PJ = F x \f\-q A *-M/* 



(49) 



and for elevation 
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PM = Fi{f\ -PJ = Wfi •«■'• A</6. (50) 

The effective quantization step size for azimuth was q A ~ 1.2 micro- 
radians, and for elevation g^ ^ 0.88 microradians, the difference being 
due to different scale factors in the two channels. The sampling period 
as stated before was At = 2 seconds. The resulting PDS of the manual 
digitizing process, P qA and P qB , are also plotted in Figs. 12 and 13. 
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Fig. 13 — Power density of random elevation angle errors, P, versus frequency, 
/. Pe , Pe\ , Pe2 — observed PDS and 95 percent confidence limits; pE.min , 
Ps.mnx = predicted tropospheric PDS limits; P t h = thermal receiver noise; 
Pqg = manual chart reading error; A; = number of degrees of freedom. 
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3.2.6 Thermal Angle Errors 

During observations of the tropospheric random angle errors it is 
important to keep angle errors due to thermal receiver noise at a com- 
paratively low level. The variance of thermal angle errors may be ob- 
tained" as 



Otk — 



K-) 



which can be reduced for S/N » 1 to 



2 

Br 



(51) 



° li = SBr(S/N) ' (52) 

where 



6 = antenna beamwidth 
S/N = input signal-to-noise power ratio 
B = receiver bandwidth 
t = post-detection integration time. 
In order to derive the power density (P lh ) of the white thermal noise 
spectrum it is first recognized that the variance of the thermal angle er- 
ror is also 

2 [ Pthdf (rn\ 

** = i 1 + (///«)' ' ( ' 

where 

f e = 1/2t (54) 

is the cut-off frequency of the post -detect ion low-pass filter. Equation 
(53) may be integrated with the substitution x = f/f c ; df = f e dx giving 

vth = /V/c-arc tan (///„) |" (55) 

or 

tr, fc 2 = |/c-P rt . (56) 

Combining (52), (54), and (56) then yields the desired density of the 
thermal angle noise as 
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0* 

Pth = 2tB(S/N) (57) 

independent of the frequency/. 

During the observations on May 7, 1965, which are analyzed in this 
report, the measured signal-to-noise ratio was (S/N)' = 23 [dB] = 
200[1] while the communications carrier of the Early Bird satellite was 
switched off. (This ratio was (S/N)' = 13 [dB] = 20 [1] due to a weaker 
beacon signal when the carrier was on.) The primed signal-to-noise 
ratios stated here are referred to a 3-kHz bandwidth. The effective noise 
bandwidth, however, is considerably lower due to the employment of a 
phase-locked tracking loop quite like the one described in Ref. 15. From 
Fig. 12 in that reference it is seen that the noise bandwidth for 

(S/N)' = 23 [dB] is B = 390 [Hz], 

which further results in an effective signal-to-noise ratio 

S/N = 200 X (3,000/390) = 1,538 [1]. 

Since the antenna beamwidth was = 0.225° = 3.94 X 10" 3 [rad] the 
desired power density of the thermal receiver noise with (57) here be- 
comes P lh = 4.1 X 10 -1 "' [rad 2 /Hz] while the communications carrier is 
switched off, and the beacon signal is strong. This thermal noise level in 
the angle measurements was low enough to permit observation of the 
random tropospheric angle variations up to frequencies of a few 0.1 [Hz], 
see Figs. 12 and 13. 

IV. COMPARISON BETWEEN PREDICTED AND OBSERVED ANGLE ERRORS 

The comparison between predicted, and observed power density spec- 
tra of random tropospheric angle variations may now be made with the 
aid of Figs. 12 and 13 into which all relevant spectra have been entered. 
The observed spectra (P A , P E ) resulted from the analysis of random 
angle error data taken on the Andover Horn to Early Bird path on May 7, 
1965 between about 23 hr:38 min EDT, and 23 hr:59 min EDT. An 
inspection of Figs. 12 and 13 shows that the PDS of the observations 
cover about two decades of frequency, namely 

2.50 X 10~ 3 [Hz] ^ / ^ 0.25 [Hz]. 

The comparison of the observed PDS (P A ; P E ) with their respective 
predicted PDS (P^.min , ^, m » s ; P*.min , Ps.max) yields almost identical 
results for the two angle oocrdinates azimuth (.4) and elevation (E). 
In particular it is found that the PDS of the observed random angle 
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variations (P A ;P E ), within their respective 95 percent confidence 
limits (P A i , Pa2 ; Pei ,Pei), He almost exactly on the predicted minimum 
power spectra (P A .min ; Pj.min) for random tropospheric angle errors. 

Thus, the observed PDS match the predicted PDS quite well in the 
shape of their frequency dependence. The low level of the observed PDS 
relative to the predicted range of PDS is thought to be due to the "good 
tracking weather" at the Andover, Maine site and at the particular time 
of observation (a quiet night). It must be remembered here that the 
prediction is based mainly upon the NBS range and phase measure- 
ments 4,5,6,7 which were obtained in Hawaii and Colorado. Whether the 
low level of the random tropospheric errors observed in Maine is a perma- 
nent property of the site, or a chance occurrence can be decided by the 
analysis of additional observations. 

It is also possible to compare the observed PDS of the azimuth errors 
with that of the elevation errors. It is found that the azimuth errors 
here have a higher level at frequencies / > 0.1 Hz than the elevation 
errors; this is an effect of the higher azimuth wind speed component 
(v A = 6.7 m/sec versus u E = 2.2 m/sec). Even larger differences between 
the azimuth and elevation random errors are expected when their ef- 
fective wind speed components differ by larger amounts. Such wind 
speed differences may be caused by either peculiar orientation of the 
natural wind vector relative to the line-of-sight, or also by differences in 
angular tracking rates. 

Near the high-frequency end of the covered band the observed PDS 
deviate significantly in shape from the predictions. This deviation is 
particularly evident in the steep increase of the observed elevation PDS 
above/ = 0.15 [Hz]. This increase is identified as an effect of quantization 
errors in the manual digitizing of the analog strip chart records. The 
transformation of these digitizing errors to the input side of the record- 
ing niters results in the steeply rising PDS (P gA , P 9 b) for these fre- 
quencies. 

In the frequency band of the observations the angle errors due to 
thermal receiver noise have a PDS (P rt ) which is negligible compared to 
that of the tropospheric angle errors, provided the communications 
carrier in the Early Bird satellite is turned off. 

It is also possible to integrate the predicted PDS of the random tropo- 
spheric angle errors over the entire error frequency band, and then to 
take the square root to obtain the standard deviation 



Pdf 
o 
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When these integrals are calculated for the predicted minimum and 
maximum PDS, it is found that the standard deviations of the random 
tropospheric angle errors are expected to lie between <r m i n PZ 10 [micro- 
radians] P& 2 [seconds of arc] and <r max PH 65 [microradians] PH 13 [sec- 
onds of arc]. This range of values compares quite well with Kennedy and 
Rosson's estimate that the tropospheric angle errors lie between 20 to 50 
microradians. 20 

The standard deviation of the expected tropospheric angle fluctuations 
versus baseline length was previously calculated from NBS data on re- 
fractivity and range variations by D. K. Barton. 21 For the equivalent 
baseline length of the Andover horn antenna of about forty feet, Barton's 
graph shows a standard deviation of perhaps seventy microradians, a 
value slightly above our predicted maximum. 

Some astronomical observations of random fluctuations in angular 
star positions, as quoted by Tatarski, 2 show standard deviations of one 
half to one second of arc. These observations have been made under con- 
ditions quite different from those for which our predictions are valid, 
namely with visible light, in clear weather, with smaller apertures, and 
probably only over a small fraction of the entire error frequency band. 
Therefore, it is not too surprising to find that these astronomical measure- 
ments lie below our minimum prediction. 

Within the limitations of the analyzed observations, and of the de- 
scribed model it is concluded that the observed random angle variations 
are essentially due to random variations of the refractive index field in 
the troposphere. The feasibility of acquiring additional data on tropo- 
spheric angle errors with the Andover horn antenna on geo-stationary 
satellites of the Early Bird type therefore was also demonstrated. These 
data may now be obtained on a routine basis with available and operat- 
ing equipment. 

The comparison of the observations given in this paper with the pre- 
diction of random tropospheric angle errors gives some confidence in the 
described analytical model. Additional observations of random tropo- 
spheric angle errors were made with radar and optical equipment over 
other propagation paths. The comparison of these observations with the 
relevant predictions from the analytical model (not reported here) are 
also satisfactory, and have further strengthened the confidence in the 
model. 

v. SUMMARY 

Earth-based radar and optical systems which are used to measure the 
position (and its time derivatives) of both distant and near objects are 
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ultimately limited in accuracy by random angle variations caused by 
fluctuations of the tropospheric refractive index. For the analysis and 
synthesis of these systems an analytical model of the random tropo- 
spheric errors has been developed. 

With this model, the predicted minimum and maximum power density 
spectra (PDS) between which observed PDS of tropospheric errors are 
expected to lie can be analytically calculated. The calculation is per- 
formed by operating with certain model functions, which depend on the 
tracking system parameters, on a model PDS (P m ) given in the range 
coordinate. P m has been derived from observations of random variations 
in the tropospheric refractive index, and in range and phase measure- 
ments made mainly at the National Bureau of Standards. 

A simplified analytical model of random tropospheric angle errors is 
described here, which is applicable to a tracking situation involving one 
(almost stationary) target and a single observer. This model is also used 
to predict the minimum and maximum PDS of the random tropospheric 
azimuth and elevation angle errors for microwave observations of the 
Early Bird geo-stationary communication satellite with the large horn- 
reflector antenna at the AT &T ground station near Andover, Maine. 

The general method of data acquisition, Fig. 5, and the specific cir- 
cumstances of some actual observations on the Early Bird satellite with 
the Andover horn are then described. Microwave azimuth and elevation 
angle measurements for an observation time of about twenty minutes 
were taken on May 7, 1965, while the Early Bird satellite appeared at 
an elevation angle of about 24.5 degrees. 

The analysis of the obtained time series of azimuth and elevation 
angles results in power density spectra (P A and P E ) and associated 
confidence limits which represent the observed random angle variations, 
see also Figs. 12 and 13. The effect of manual chart reading errors on the 
observed PDS was also studied. It was shown to consist of a steep in- 
crease in the PDS at the high frequency end. The effect of thermal 
receiver noise on the observed random angle variations was kept at a 
negligible level. 

The comparison of the predicted PDS of the random tropospheric 
angle errors for the Early Bird observations with the observed PDS leads 
to the conclusion that the observed random angle variations are indeed 
caused by the troposphere. In particular it is found that the PDS of the 
azimuth and elevation observations (P A and P E in Figs. 12 and 13), 
within their respective confidence limits (P A i ,Pa2 ; Ps\ , P«), lie almost 
exactly on the predicted minimum power density spectra (P A , min ; 
Pr, min) for random tropospheric angle errors. 
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The feasibility of acquiring additional data on tropospheric propaga- 
tion effects, especially random angle errors, with the Andover horn 
antenna on geo-stationary satellites of the Early Bird type, therefore, 
was also demonstrated. 
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APPENDIX 

Numerical Calculation of Predicted PDS 

The parameters of the tracking situation during the Early Bird observa- 
tions on May 7, 1965, which permit the prediction of the tropospheric 
angle PDS with the described model are: 

transmission frequency = ft = 4137.86 [MHz] 

antenna diameter = d = 67.7 [ft] = 20.6 [m] 

beamwidth = = 0.225 [deg] 

apparent elevation angle = E a re 24.5 [deg] 

azimuth angle = A re 128.5 [deg] 

altitude of horn antenna = fa = 900 [ft] = 274 [m] 

altitude of satellite = fa re 22,200 [st. mi.] re 35,700 [km] 

slant-range = Rn re 24,300 [st. mi.] re 39,100 [km] 

wind vector: | 4> \ = 19 [st. mi/hr]; 5 = 0° 

surface refractivity = N, = 301 

With these parameters the model functions for this tracking situation 
are calculated as follows. 
Breakfrequency of the angle scale function, (15): 

/i = 2.43 X 10~ 2 [Hz]. 

Angle scale function, (14): 

S a = 20 (//Hz) 2 (l/m 2 ] for ^ / ^ 2.43 X 10 -2 [Hz] 
S a = 1.18 X 10~ 2 [1/m 2 ] for 2.43 X 10~ 2 [Hz] ^ / ^ oo. 

Breakfrequency for aperture smoothing, (17): 
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f 2 = 9.71 X 10 -2 [Hz]. 

Aperture smoothing function, (16) : 

$ a = 1 for ^ / £ 9.71 X 10" 2 [Hz] 

<*>« = 9.43 X 10~ 3 (//Hz) -2 for 9.71 X 10" 2 [Hz] £ / £ oo. 

Effective tropospheric path length, (19) and (20): 

L = 8.02 [km]. 

Effective path length function, (18) : 

A = 0.534 for = / ^ « . 

Weather functions: 

minimum: W mi „ = 1 (by definition of P m ) 

maximum: TF ma x as per (22). 

Effective wind speed, azimuth, (30): 

u A = 14.9 [st. mi/hr] = 6.7 [m/sec]. 

Effective wind speed, elevation, (33 ) : 

u E = 4.9 [st. mi/hr] = 2.2 [m/sec]. 

Effective wind function for transformation of the power density, azi- 
muth, (26): 

U PA = 0.149 7+1 . 

Effective wind function for transformation of the frequency regions, 
azimuth, (27): 

U tA = 6.7. 

Effective wind function for transformation of the power density, ele- 
vation, (26): 

Un = 0.455 T+1 . 

Effective wind function for transformation of the frequency regions, 
elevation, (27): 

U/m = 2.2. 

The operation with these model functions on the range model PDS 
P m \j) results in the following four predicted PDS of Iropospheric random 
angle errors. 
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Minimum PDS in azimuth: 

(7.57 X 10 +16 (//Hz)"* [rad'/Hi] 

for ^ / ^ 1.68 X 10~ 7 [Hz] 

3.55 X 10 -4 (f/Hz) +1 [rad 2 /Hz] 

for 1.68 X 10~ 7 ^ / ^ 6.70 X 10 6 [Hz] 

1.95 X 10 -10 (//Hz) -05 [rad 2 /Hz] 

for 6.70 X 10" B ^ / ^ 6.70 X 10 3 [Hz] 

1.60 X 10"" (//Hz) -1 [rad 2 /Hz] 
for 6.70 X 10 -3 ^ / ^ 1.63 X 10 _1 [Hz] 

4.25 X 10" 13 (//Hz)" 3 [rad 2 /Hz] 

for 1.63 X 10" 1 £ / £ 6.51 X 10" 1 [Hz] 

1.81 X 10~ 13 (f/Hz) -5 [rad 2 /Hz] 

for 6.51 X 10" 1 £ / £ 6.70 X 10 +2 [Hz] 

5.47 X 10 -6 (//Hz) -8 [rad 2 /Hz] 

for 6.70 X 10 +2 £ / ^ * [Hz] 

Maximum PDS in azimuth: 

f4.54 X 10 +17 (//Hz) +4 [rad 2 /Hz] 

for ^ / g 1.68 X 10" 7 [Hz] 

2.13 X 10" 3 (f/Hz) +1 [rad 2 /Hz] 

for 1.68 X 10 -7 ^ / ^ 1.49 X 10 4 [Hz] 

3.90 X 10" 9 (f/Hz) -0 " 5 [rad 2 /Hz] 

for 1.49 X 10" 4 £ / £ 1.63 X 10" 1 [Hz] 

1.03 X 10- 10 (//Hz)" 2 6 [rad 2 /Hz] 

for 1.63 X 10 _1 g / ^ 6.51 X 10 _1 [Hz] 

4.38 X lO -11 (//Hz) -4 5 [rad 2 /Hz] 

for 6.51 X 10 -1 ^ / ^ 6.70 X 10 ' [Hz] 

3.61 X 10"" (//Hz)" 5 [rad 2 /Hz] 
for 6.70 X 10 _1 ^ / ^ 6.70 X 10 +2 [Hz] 

1.09 X 10" 2 (//Hz) -8 [rad 2 /Hz] 

for 6.70 X 10 +2 ^ / ^ co [Hz] 



A ,max 



" B.min — 
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Minimum PDS in elevation: 

2.01 X 10 +19 (//Hz) +4 [rad 2 /Hz] 

for ^ / ^ 5.50 X 10" 8 [Hz] 

3.31 X 10" 3 (//Hz) +1 [rad 2 /Hz] 

for 5.50 X 10~ 8 ^ / ^ 2.20 X 10" 6 [Hz] 

3.41 X 10" 10 (/"/Hz) -0 - 8 [rad 2 /Hz] 

for 2.20 X 10~ 5 g / g 2.20 X 10~ 3 [Hz] 

1.60 X 10~ u (//Hz)" 1 [rad 2 /Hz] 

for 2.20 X 10" 3 ^ / ^ 5.35 X 10" 2 [Hz] 

4.57 X 10~ 14 (//Hz)" 3 [rad 2 /Hz] 

for 5.35 X 10" 2 g / ^ 2.14 X 10 _I [Hz] 

2.09 X 10" 15 (//Hz)" 5 [rad 2 /Hz] 

for 2.14 X 10 _1 g / g 2.20 X 10 +2 [Hz] 

2.22 X 10~ 8 (//Hz)" 8 [rad 2 /Hz] 

for 2.20 X 10 +2 ^ / ^ oo [Hz] 

Maximum PDS in elevation: 

fl.21 X 10 +2 ° (//Hz) +4 [rad 2 /Hz] 

for ^ / ^ 5.50 X 10" 8 [Hz] 

1.99 X 10" 2 (//Hz) +I [rad 2 /Hz] 

for 5.50 X 10" 8 g / ^ 4.91 X 10" 5 [Hz] 

6.81 X 10" 9 (//Hz)" 05 [rad 2 /Hz] 

for 4.91 X 10" 6 g / g 5.35 X 10" 2 [Hz] 

1.95 X 10" 11 a/Hz) -25 [rad 2 /Hz] 

for 5.35 X 10~ 2 ^ / ^ 2.14 X 10" 1 [Hz] 

8.91 X 10" 13 (//Hz)" 4 5 [rad 2 /Hz] 

for 2.14 X 10 _1 g / ^ 2.20 X 10 -1 [Hz] 

4.17 X 10" 13 (//Hz)" 5 [rad 2 /Hz] 

for 2.20 X 10" 1 g / ^ 2.20 X 10 +2 [Hz] 

4.43 X 10" 6 (//Hz)" 8 [rad 2 /Hz] 



*E,nt 



for 2.20 X 10 + " g / g oo [Hz] 
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